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The BCS-BEC crossover within the two-dimensional attractive Hubbard model is studied by using 
the Cellular Dynamical Mean-Field Theory both in the normal and superconducting ground states. 
Short-range spatial correlations incorporated in this theory remove the normal-state quasiparticle 
peak and the first-order transition found in the Dynamical Mean-Field Theory, rendering the normal 
state crossover smooth. For U smaller than the bandwidth, pairing is driven by the potential energy, 
while in the opposite case it is driven by the kinetic energy, resembling a recent optical conductivity 
experiment in cuprates. Phase coherence leads to the appearance of a collective Bogoliubov mode 
in the density-density correlation function and to the sharpening of the spectral function. 
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The problem of the crossover between the BCS and 
the Bose-Einstein condensation (BEC) has been of great 
interest in the context of the pseudogap observed in un- 
derdoped cuprates [|. The recent discovery of the BCS- 
BEC crossover in ultracold fermionic atoms trapped in 
optical lattices |2| has renewed our interest in this is- 
sue. Does the condensation of loosely bound Cooper 
pairs evolve smoothly into the BEC of tightly bound 
composite bosons as the attraction between fermions is 
increased gradually? This question was first addressed 
by Leggett Q who proposed that these two pictures are 
limiting cases of a more general theory in which both the 
fermionic nature of individual particles and the bosonic 
nature of pairs must be considered on an equal footing. 
By using a T— matrix approximation in the intermediate 
coupling regime, P. Nozieres and Schmitt-Rink Q ex- 
tended Leggett 's analysis to a lattice and to finite tem- 
perature to find a smooth connection between the two 
limits. Since then, various other approximate schemes 
have been used to understand the pseudogap phenomena 
as well as the BCS-BEC crossover. For example, Quan- 
tum Monte Carlo (QMC) simulations HH0 of the two- 
and three-dimensional attractive Hubbard model found 
a breakdown of Fermi liquid theory in the normal state, 
accompanied by a pseudogap and a spin gap. 

The first Dynamical Mean-Field Theory (DMFT) 
study of the BCS-BEC crossover at arbitrary attractive 
interactions was carried out by Keller et al. @- The 
authors calculated the transition temperature for super- 
conductivity, which smoothly interpolates from the BCS 
behavior at weak coupling to the t 2 /U behavior at strong 
coupling. However, the double occupancy, the spin sus- 
ceptibility, and the quasiparticle weight show, in the 
T — > limit of their normal state solution, a disconti- 
nuity near U ~ 1.5W with W the bandwidth. Capone et 
al. 10] studied the first-order transition in detail by using 
exact diagonalization as the impurity solver for DMFT 
and later extended their work to finite temperature [lT| . 
The discontinuity found in the normal state of DMFT 



may suggest a radically different mechanism for super- 
conductivity at weak and strong coupling, but that is 
problematic in the context of the expected smooth BCS- 
BEC crossover (in the superconducting state). Neverthe- 
less, in a recent optical conductivity experiment in the 
cuprates, Deutscher et al. [T^ found that near optimal 
doping there is a reversal of the sign of the kinetic energy 
difference between the superconducting (SC) and normal 
(NR) states. Although the mechanism for superconduc- 
tivity in the cuprates is still under debate, this intriguing 
experiment calls for explanation also in the context of 
the BCS-BEC crossover, as the authors noted. 

In the BCS limit, it is well established that supercon- 
ductivity is potential-energy driven since the broadening 
of the Fermi surface caused by superconductivity leads 
to an increase in kinetic energy in the superconducting 
state. At strong coupling, in the BEC limit, one may 
expect that superconductivity is kinetic- energy driven 
based on two different arguments : a) The gain in poten- 
tial energy occurs at high temperature where the bosons 
form. At low temperature, to leading order in a low den- 
sity expansion [F| , the bosons condense for the same rea- 
son as free bosons, because of the gain in kinetic energy. 
Hence, comparing with a normal state where the bosons 
are formed, the superconductivity occurs because of gain 
in kinetic energy, b) One can map the attractive Hubbard 
model to the half-filled repulsive model in the presence of 
a magnetic field. The antifcrromagnetic order that ap- 
pears close to half-filling in this model is analogous to 
superconductivity and at strong coupling, it occurs be- 
cause having nei ghbo rs that are ordered leads to a gain in 
exchange energy hj. One can check that in the mapping 
of the Hubbard model to the Heisenberg or t-J model, the 
exchange energy corresponds to minus twice the poten- 
tial energy of the original Hubbard model These 
physical arguments do not tell us whether the change in 
pairing mechanism occurs in a continuous or discontin- 
uous manner, at what coupling they occur or what is 
the order of magnitude of the condensation energy. The 
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explicit calculations presented in the present paper thus 
answer both qualitative and quantitative questions. 

First we emphasize our main results: 1) Including 
short-range spatial correlations explicitly removes the 
first-order transition found in DMFT, rendering the 
crossover in the NR state also smooth, without a dis- 
continuity, just as in the SC state. 2) Near U equal to 
the bandwidth of 8i, a change in pairing mechanism oc- 
curs. For U < 8t the condensation energy is lowered by 
the potential energy while for U > 8t it is lowered by the 
kinetic ener gy, resembling a recent optical conductivity 
experiment |l2j. 3) The phase coherence manifests itself 
most dramatically by the appearance of a collective Bo- 
goliubov mode in the density-density correlation function 
and by the sharpening of the spectral function compared 
with the normal state. 

Using Cellular Dynamical Mean-Field Theory 
(CDMFT) 0, we study the crossover between weak 
and strong-coupling in the two-dimensional attractive 
Hubbard model [Hf 

where c ia (ci a ) are creation (annihilation) operators for 
electrons of spin a, rn a — <y ia Ci a is the density of a spin 
electrons, t$j is the hopping amplitude equal to — t for 
nearest neighbors only, — U is the on-site attractive in- 
teraction with U > and /i is the chemical potential 
controlling the electron density. The CDMFT method 
is a natural generalization of the single site DMFT 8] 
that incorporates short-range spatial correlations. In the 
CDMFT construction [lii llq the infinite lattice is tiled 
with identical clusters of size N c , and the degrees of free- 
dom in the cluster are treated exactly while the remaining 
ones are replaced by a bath of non-interacting electrons 
that is determined self-consistently. Since the CDMFT 
method treats short-range spatial correlations explicitly, 
it is able to describe features caused by finite dimension- 
ality or finite coordination number, which are missed in 
the single site DMFT. 

To solve the quantum cluster embedded in an effective 
SC medium, we consider a cluster-bath Hamiltonian of 
the form [12, HJ 

"f" ^ ^ £ m<T a rn<r a ma- ~F ^ , ^rap,u ("LffV "F H.C.) 

+ ^A(a« T a« i+ H.c.). (2) 

Here the indices fi, v = 1, • • • ,N C label sites within the 
cluster, and and a^ na annihilate electrons on the clus- 
ter and the bath, respectively. In the present study we 
used N c = 4 sites for the cluster (minimum number 
of sites reflecting the full square lattice symmetry) and 
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FIG. 1: (a) Uniform static cluster spin susceptibility \sp and 
double occupancy rid =< n^riii >. (b) Imaginary part of the 
s— wave cluster pair correlation function with total momen- 
tum equal to zero Xsc(w). The above quantities are calculated 
in the NR ground state at quarter filling (n = 1/2). 



Nb = 8 sites for the bath with m = 1, 4, a = 1, 2. t^v 
is the hopping matrix within the cluster and, using sym- 
metry, e^ a = e a is the bath energy and = V a S m ^ 
is the bath-cluster hybridization matrix. A represents 
the amplitude of s— wave SC correlations in the bath. To 
deal with superconductivity, the Nambu spinor represen- 
tation is used for the cluster operators 

* t = ( c It> C 2T> C 3T' C 4T> C H> C 217 C 31,C4x) , 

while the Weiss field, the cluster Green's function and 
self-energy constructed from these operators are 8x8 ma- 
trices. The exact diagonalization method [2lJ is used to 
solve the cluster-bath Hamiltonian Eq.[3at zero temper- 
ature, which has the advantages of computing dynamical 
quantities directly in real frequency and of treating the 
large U regime without difficulty. All the results pre- 
sented here are obtained at quarter filling (n = 1/2) 
far from the particle-hole symmetric case where charge- 
density- wave and pairing instabilities coexist. Qualita- 
tively similar results were found at other densities. 

Fig.^a) presents the evolution with U of the uniform 
static cluster spin susceptibility \sp and of the double oc- 
cupancy rid in the NR state obtained by forcing SC order 
to vanish. In the case of the attractive Hubbard model, 
Xsp and rid describe how many fermions turn into local 
singlet pairs due to attraction. In the limit of U — » oo, 
Xsp — ► and rid —> n/2, while in the opposite limit, 
Xsp -> 2N(0) and n d -> {n/2) 2 . Here N(0) is the (non- 
interacting) density of states per spin. ^From weak to 
strong coupling these thermodynamic quantities are con- 
tinuous in CDMFT, in stark contrast to those of the sin- 
gle site DMFT (see Fig.2 and Fig.3 in Ref. @). Near 
U/t = 8 both 

Xsp an d rid start to saturate, indicating 
that tightly bound bosonic pairs 01 begin to dominate 
the physics. 

Fig. |3 shows additional evidence of the absence of a 
first order transition when short-range correlations are 
treated explicitly. Unlike in DMFT, where the spectral 
function A(k,ui) has a peak at uj = that disappears at 
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FIG. 3: (a) Kinetic Ek, potential Eu and total Et energies 
computed in the SC and NR ground states denoted as the 
solid and dashed curves, respectively, (b) Total energy differ- 
ence between the SC and NR states AEt- The condensation 
energy obtained from the BCS theory is shown as diamonds. 
The circles represent T c obtained by Keller et al. |i| within 
DMFT but multiplied by 0.09. 



FIG. 2: Single particle spectral function A[k,ui) at k = 
(3/87T, 3/87r) near the Fermi surface calculated in the NR 
ground state with a broadening parameter of 0.125t. The 
dashed curve in the first panel is computed from the TPSC 
at finite temperature of (3 = t/T = 8 on a 64 x 64 cluster. 



a critical coupling U c /t ~ 12 [23 to lead to an insula- 
tor, in CDMFT A(k,u>) is a minimum at w — already 
at U/t = 6 and we claim that the apparent absence of 
a gap at U/t = 4 (the first panel) is an artifact of the 
finite size of the cluster used. In the case of weak cou- 
pling, short-range correlations available in a cluster of 
size N c = 4 are not long enough to lead to a gap, as 
was found in the repulsive Hubbard model |2^|. When 
a large 64 x 64 lattice is used to capture long-range cor- 
relation effects within the Two-Particle Self-Consistent 
(TPSC) theory [24| (dashed curve) valid at weak cou- 
pling, a gaplike feature does exist even at a finite tem- 
perature of T/t =1/8. QMC calculations on a 16 x 16 
lattice show the same result. Hence, there is enough 
evidence that as soon as finite-range spatial correlations 
are included explicitly, the first order transition from a 
Fermi liquid to a non-Fermi liquid state found in the NR 
state solution of the DMFT equations disappears imme- 
diately [2f| , making the NR state crossover smooth [2(| , 
just as in the SC state. Physically, the reason why DMFT 
leads to a first-order transition is because residual hop- 
ping between preformed pairs is of order t 2 /U xl/d with 
d the dimension, hence they have vanishing kinetic en- 
ergy contribution in large dimension and they localize. 
Within CDMFT boson hopping is restored (within the 
clusters). With increasing U, A(k,u>) has a large in- 
coherent spectrum at high frequencies together with a 
sharp peak near the gap edge. A close inspection shows 
a continuous evolution of high energy peaks for U/t > 6. 

The smooth crossover can also be seen by consider- 
ing Fig. ^b), which shows the imaginary part of the 



s— wave cluster pair correlation function computed in the 
NR state for several couplings. With increasing U be- 
yond 6t the peak intensity increases, while the peak po- 
sition decreases, scaling as J. This is analogous to the 
t 2 /U scaling of the peak position in the spin-spin corre- 
lation function at q = (ir, it) in the half-filled Hubbard 
model 0| . The gradual change of the pair correlation 
function from weak to strong coupling is consistent with 
the smooth crossover of A(k,u) as a function of U. 

Next we study the BCS-BEC crossover in the SC state 
and its consequences. Fig. [31(a) shows the kinetic Ek, 
potential Eu and total Et energies computed in the SC 
and NR ground states denoted as the solid and dashed 
curves, respectively. Generally the energy differences are 
tiny (only about 1 — 2%), which implies that the two 
ground states are energetically very close. At weak cou- 
pling one finds, as in the BCS theory [2^|, that while 
kinetic energy is increased in the SC state by the broad- 
ening of the Fermi surface, the decrease in the potential 
energy overcompensates it. At strong coupling, however, 
the roles are interchanged. On a lattice, phase coher- 
ence involves breaking up of more local pairs by virtual 
hopping (increase in the potential energy) to enhance 
their mobility further (decrease in the kinetic energy). 
The full inversion of the roles occurs near U/t = 8, but 
the decrease in kinetic energy happens somewhat earlier 
(U/t = 6) in our study. Although the mechanism for 
high temperature superconductivity is debated and the 
observed gap symmetry (d— wave) is different from the 
s— wave predicted in the attractive Hubbard model, the 
current result resembles a recent optical conductivity ex- 
periment in the cuprates 01 where AEk crosses over 
from a BCS behavior (AEk > 0) to an unconventional 
behavior (AEk < 0) as the doping decreases. In both 
the SC and the NR states, the kinetic energy scales as 
t 2 /U. The difference in the total energy AEt is plotted 
in Fig. Ofb) as a solid curve. The condensation energy 
scales with the T c found by Keller et al. 9] in DMFT 
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FIG. 4: (a) Imaginary part of the collective cluster charge 
excitations Xchiui) at q = in the SC and NR (inset) ground 
states. The solid, dotted, dashed and long-dashed curves cor- 
respond to U/t =4, 6, 8 and 12, respectively. Xch{u>) for 
U/t —S and 12 in the inset are invisibly small on this scale, 
(b) Density of states N(u>) computed in the SC (solid) and 
NR (dashed) states at U/t = 8. A broadening parameter of 
0.125£ is used in (a) and (b). 

(circles), namely it reaches a maximum near U/t = 8 
and decreases as t 2 /U beyond it. At U/t = 4 it is al- 
ready far smaller than that in the BCS theory. Despite 
the identical t 2 /U scaling of AEt in our result and of 
T c in DMFT at strong coupling, they differ by a factor 
of 10 approximatively. This is a upper bound since one 
expects the mean-field like DMFT estimate for T c to be 
high because of the neglect of spatial fluctuations. 

As was shown above, phase coherence does not lead 
to a large gain in energy. The difference between the 
NR and SC states should manifest itself most dramati- 
cally in correlation functions. We computed the uniform 
charge-charge correlation function in the cluster in both 
states, as shown in Fig. EJa). With increasing U the 
mode energy decreases, in sharp contrast with the NR 
state spectrum shown in the inset. For U/t > 6 in the 
SC state, the charge excitation forms a sharp resonance 
much below twice the gap, while in the NR state (inset) 
the weight moves away from u> — with weight around 
lu = U (not shown) that corresponds to breaking a pair. 
In fact it is a collective Bogoliubov mode [2^| which is sep- 
arated from the continuum that comes from the breaking 
of Cooper pairs at higher energy (not shown) . Because a 
local pair becomes well-defined only at strong coupling, 
the difference between the two states is most pronounced 
for U/t — 8 and beyond. The difference between the two 
states also appears in the local density of state N(u) but 
is more subtle, as shown in Fig.^Jb). Although a pairing 
gap without phase coherence already exists in the NR 
state, the condensation depletes the remnant spectrum 
from the gap and builds a sharp spectrum at the gap 
edge. 

Note that in our approach the phase fluctuations are 
short range. Hence, long-wavelength fluctuations do not 
appear necessary to lead to kinetic-energy driven pairing. 
Also, CDMFT allows one to go beyond Eliashberg the- 
ory. The latter is a "strong coupling" theory in a different 



sense than that used in the present paper. It cannot treat 
the case where U is of order the bandwidth or larger, 
i.e. the case where the normal state consists of bound 
fcrmions. The present study focuses on a 2 x 2 cluster 
with 8 bath sites since the next smallest cluster size that 
corresponds to 0.5 filling would need unreasonable com- 
puter resources. Nevertheless low (but finite) tempera- 
ture CDMFT+QMC calculations lead us to believe 
in the robustness of our results except at the weakest 
coupling. 

The increasing fractional deviation of the DMFT [3l| 
and the CDMFT order parameter and the energy gap 
from their corresponding BCS values with decreasing U 
may be caused not only by the absence of quantum fluc- 
tuations in the BCS theory, but also by the fact that 
clusters (in DMFT and CDMFT) are smaller than the 
spatial extent of the Cooper pairs in that limit. 

Although short-range spatial correlations incorporated 
in this theory remove the first-order transition found in 
the NR state of DMFT, large frustration may lead back 
to a first-order transition even in CDMFT as was found in 
the half-filled repulsive Hubbard model. For that model, 
we found for t'/t = -0.717 (t'/t = 1 in the BEDT model) 
that a first-order transition does occur at U close to the 
bandwidth 27]. In the attractive Hubbard model, frus- 
tration could come from magnetic field for example. 

To summarize, we have studied the BCS-BEC 
crossover within the two-dimensional attractive Hubbard 
model by using the Cellular Dynamical Mean-Field The- 
ory (CDMFT) both in the normal and superconducting 
ground states. Explicit treatment of short-range spatial 
correlations in this theory removes a first-order transi- 
tion found in DMFT, making the normal state crossover 
smooth without a discontinuity. For U smaller than the 
bandwidth, pairing is driven by the potential energy, 
while in the opposite case it is driven by the kinetic en- 
ergy, resembling a recent optical conductivity experiment 
in the cuprates. The condensation energy has a maxi- 
mum at U equal to the bandwidth and scales as t 2 /U 
at strong coupling. Phase coherence leads to the ap- 
pearance of a collective Bogoliubov mode in the density- 
density correlation function and to the sharpening of the 
gap already present in the normal state spectral function. 

Note added in proof: Recently we became aware of 
DMFT results by Toschi, Capone, and Castellani [32| . 
which also showed potential-energy to kinetic-energy- 
driven pairing with increasing U in the single site DMFT. 
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